Fermionic concurrence in the extended Hubbard dimer 
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In this paper, we introduce and study the fermionic concurrence in a two-site extended Hubbard 
model. Its behaviors both at the ground state and finite temperatures as function of Coulomb 
interaction U (on-site) and V (nearest-neighbor) are obtained analytically and numerically. We also 
investigate the change of the concurrence under a nonuniform field, including local potential and 
magnetic field, and find that the concurrence can be modulated by these fields. 
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I. INTRODUCTION 



Recently, many efforts have been devoted to the entan- 



trivial behavior in these systems may shed new hghts on 
the understanding of physical phenomena of condensed 
matter physics. A typical case is the relation of entangle- 
ment to quantum phase transitionjQ, Q, IM 1^ Hy; HB 113 ■ 
For example, Osterloh el a/., pi reported that the entan- 
glement shows scaling behavior in the vicinity of quan- 
tum phase transition point of the transverse-field Ising 
model. Most of previous works are for spin 1/2 systems, 
where the degrees of freedom at each site (or qubit) is 2. 
For these sy sterns, the entanglement of formation, i.e., 
the concurrence |l7l| . is often used as a measure of pair- 
wise entanglement between two qubits. Such measure is 
only valid for 2 (g) 2 systems. If the degrees of freedom at 
each qubit is larger than 2 (for example, the spin 1 sys- 
tem or systems consisting of fermions with spin), how to 
quantity the entanglement of arbitrary entangled state is 
a challenging issue. Several studies ,18j, 19, 20, 21] were 
devoted to this issue. For example, Mintert et al. ob- 
tained a lower bound for the concurrence of mixed bipar- 
tite quantum states in arbitrary dimensions. Neverthe- 
less, it is still very difficult to provide a reliable measure 
for the pairwise entanglement of systems with the number 
of local states larger than 2. To the best of our knowl- 
edge, none of previous work investigated the pairwise en- 
tanglement for systems consisting of electrons with spin, 
such as the Hubbard model, although there were a few 
works studied the local entanglement of fermionic models 
Ll4, 15, 16]. 

In this paper, we introduce and study the fermionic 
concurrence by using the extended Hubbard dimer as 
an example. Besides its importance in exploring many- 
body correlation in condensed matter physics, a dimer 
system also has potential utility in the design of quantum 
device 22] . By considering the Pauli's principle, we first 
illustrate how to quantify the fermionic concurrence in 
the Hubbard model and formulate it in terms of fermion 
occupation and correlation so one can easily calculate it. 
Then based on the exact solution of the Hubbard dimer, 
we obtain the result at the ground state and show that 



the fermionic concurrence could be used to distinguish 
state exhibiting charge-density correlation from state ex- 
hibiting spin-density correlation. We also study its be- 
havior at finite temperatures |2] , and find that it is a de- 
creasing function of temperature. Moreover, we investi- 
gate the behavior of the concurrence under a nonuniform 
local potential and magnetic field 23, 24|[23- We find 
that the concurrence could be modulated by these local 
fields. Our work therefore not only provides a possible 
way to investigate the pairwise entanglement in the elec- 
tronic system, but also enriches our physical intuition on 
the systems with both charge and spin degree of freedom. 
Some results are instructive for future experiments. 



II. THE MODEL AND FORMULISM 

The Hamiltonian of the one-dimensional extended 
Hubbard model reads 

H = -tY^ c]^^Cj+s,a + U^nj^nji + V^njnj+i.{l) 

where a =t, i; j ^ 1, . . . ,L;6 — ±1, c'- and Cja create 
and annihilate an electron of spin a at site j, respec- 
tively, and the hoping amplitude t is set to unit. At each 
site, there are four possible states, \0)j, \ t)j, | i)j, \ T-l)j 



denoted by 



/]' 



V — 1,2,3,4. The Hilbert space of 



L-site system is of 4 dimension, and |i^i, V2 ■ ■ -v^) 



n 



are its natural bases. Therefore any state in 



such a system can be expressed as a superposition of 
the above bases. We consider reduced density matrix 
Pji — Tijip of site j and I, where p is the thermal density 
matrix of the system, and Tr^-; stands for tracing over 
all sites except the freedom of jth and Ith site. Thus 
Pji defines the quantum correlation between site j and 
site /. However, since there does not exist a well-defined 
entanglement measurement for a mixed state of 4 (g) 4 
bipartite systems, it is impossible to study the entangle- 
ment between two sites exactly. Fortunately, the Hilbert 
space of site j and I can be expressed in a direct-product 
form of electrons with spin up and spin down, that is, 
for site j and /, we have two subspaces, one with the 
bases |0,0), |0, t), | T^O),] tiT)i and the other one with 
|0, 0), |0, i), I i, 0), I i, i), respectively. It is therefore pos- 



sible to investigate the problem in a given sector with 
spin up or spin down separately. For convenience, we re- 
strict our studies in the subspace of spin up electrons, and 
the other part can be simply obtained via the Z2 symme- 
try. The Hamiltonian |^ possesses U(l)xSU(2) symme- 
try, i.e., it is invariant under the gauge transformation 
Cj.a —> e^^Cj,a and spin rotation Cj^a —^ UaSCj^s, which 
manifest the charge and spin conservation. The latter 
implies the absence of coherent superposition between 
configurations with different eigenvalues of tij^^ -\- rij^i-^. 
Thus the reduced density of electrons with spin-up on 
two sites has the form 



Pj,i 



fu'^ Q 

wi z 

Z* W2 

u- 



(2) 



in the standard bases |0, 0), |0, t), | T: 0)j I T: T)- Elements 
in the density matrix pji are related to single occupation 
and correlations between the two sites. 



,]CIA) 






u 



l-2(nj,|) + (nj,|nj+i,|). 



(3) 



where () denotes the expectation value of the correspond- 
ing operator. 

We use the concurrence as a measure of entangle- 
ment for such two-qubit system. It is defined in terms 
of the spectrum of the matrix Pj/Pj/jlJl where pji 



V 

crfp 



PjiPji and Ai > A2 > A3 > A4 
be calculated as 



C = max 0, v '^i 



Precisely, if A^s are eigenvalues of 
the concurrence can then 



(4) 



Since there exists a monotonous relation between the con- 
currence and the entanglement of formation Ef^ Ef = 
— a:log2a;— (1 — a;)log2(l — x), where x ~ l/2 + \/l — C^/2 
|l7l | , we will hereafter use the concurrence instead of en- 
tanglement of formation in our study. From Eq. (|2Jl, the 
fermionic concurrence can be calculated as 



C ~ 2 max 



{k4t 



Ci,T/ 



^i.T";. 



), 0} 



(5) 



III. HUBBARD DIMER WITH TWO 
ELECTRONS 

In this section, we consider a model which consists of 
two sites and two electrons, because not only can it be ex- 
actly solved, but also it gives us a clear physical picture. 
The Hamiltonian for the dimer reads 

(7 

-|-C/(ni,|ni,x -t- n2,-\n2,i) + 2Vnin2. (6) 



In the standard bases | T,i),| i,T),l Ti,0),|0,Ti) of the 
reduced subspace with zero magnetization: S^ = Q, H is 
a 4 X 4 matrix 



H^ 




(7) 



and we can easily obtain its eigenfunction 

*i = ^ [i n,o) + |o,n)] + ^ [| ta) + 1 u)] 

vi/2 = ^ [i n,o) + |o,n)] - ^ [| T,i> + 1 i,T)] 
1 



*3 = ^[in,o)-|o,n) 

*4 = ^[|T,i>-U,T)] 



(8) 



where 9 is defined as 
tan( 



(C//2 -V) + yJ{U/2 - Vf + 16 
The corresponding eigenvalues are 

El = ]^{U + 2V- ^M+{U-2VY), 
1 



E2 = -{U + 2V + v/64+(f/-2T/)2), 

Ei ^ 2V. 
^From the ground-state wavefunction, we have 



(9) 



{^i tC2,y — smt^cost^ 
(ni,T"2,T> = 0. (10) 

Thus, the fermionic concurrence takes the form of 



C = |sin26l|. 



(11) 



Clearly, the ground-state wavefunction simultaneously 
comprises the states of the double occupied state | t| 
, 0) -I- |0, ti) and the single occupied state | T, i) + I ■!,?)• 
The magnitude of d determines which of them dominates 
in the ground state. From its definition, it is easy to 
learn that the line U = 2F, i.e. — 45°, separates the 
two regions. One region dominated by | tijO) -I- |0, ti) 
corresponds to the charge-density-wave state for a large 
system, and the other one dominated by | Tii) + I i 
, t) corresponds to the spin-density-wave state. Clearly, 
the fermionic concurrence reaches a maximum at the line 
V — C//2, and shows the mirror symmetry along the line 
y = C//2 (See figure UJ. 

At finite temperatures, the expectation value of a 
given operator O is (O) = ^n^n\0\n)e~f^^'^ jZ, where 
/? = 1/T, Z is the partition function, and E^ denotes the 





FIG. 1: The fermionic concurrence of the extended Hubbard 
dimer as a function of U and V. The curves in the C — 
plane constitutes a contour map, and the thick hne denotes 
the local extremum. 




FIG. 2: The fermionic concurrence of the Hubbard dimer 
{V = 0) as a function of U and T. 



energy spectra. For the present model, we should also 
take into account states in the subspace with S*^ = ±1 
besides that of S^ — 0. Then we have two additional 
states ^5(6) = I TiDd iji)) with the same eigenvalue 
EF.tf^\ = 2V. Thus the fermionic concurrence of the Hub- 



5(6) 

bard dimer takes the form. 



C(T) = -max[|sin26l(e-''^i 



-liE^i 



)|-2e- 



-!iE^ 



0](12) 



It is not difficult to prove that the above expression is 
a decreasing function of T, as shown in figure [3 for a 
specific value of T^ = 0. 

In the high temperature limit, T —f oo, the Boltzmann 
weight of each eigenstate becomes almost equal, so the 



FIG. 3: The threshold temperature Tth of the fermionic con- 
currence as a function of U and V. 



first term in the square brackets of Eq. (|12|l becomes 
negative and we have a vanishing fermionic concurrence. 
Hence we expect that there exists a threshold temper- 
ature Tth at which the concurrence becomes zero. Pre- 
cisely, it is determined by 



sin 20{e 



-Bi/Tth 



--Ea/Tth^ 



) - 2e"-^=/^"^ = 0. (13) 



Clearly, the hoping-correlation (cji^c/j) arises from the 
state of ^1 and ^2, and charge correlation (ni,|n2,t) from 
^5. Since the Boltzmann weight e"^^/^"* of ^1 is always 
larger than e"^^/"'""', the term of absolute value is then 
positive for finite T. However, the term U + 2V in the 
expression of Boltzmann wight breaks the mirror sym- 
metry of the fermionic concurrence at the ground-state 
along the line U — 2V. For example, the concurrence in 
the fourth quadrant of the U — V plane is more susceptive 
to thermal fluctuation and might be easily suppressed to 
zero than other quadrants. Therefore, in different quad- 
rants, the effect of thermal fluctuation is quite different. 
Then from Eq. (|13|l . we obtain a slope- like threshold 
temperature for the concurrence (see figure |3J). 



IV. UNDER A NONUNIFORM FIELD 

In this section, we would like to consider the behavior 
of the concurrence under an external field. In general, 
a uniform external field, like the local potential or the 
magnetic field, suppresses the concurrence since it frus- 
trates the hoping process. Thus we are only interested 
in the case with a nonuniform field. Moreover, differing 
from the previous investigation where calculations were 
done for canonical ensembles, we perform calculation for 
grand canonical ensembles. The first case we consider is 




chemical potential: u 



FIG. 4: The fermionic concurrence as a function of the local 
potential /x (u instead of ^ as an indictor for axes) and V at 
the ground-state. Here, U = 2. 




FIG. 6: The fermionic concurrence as a function of magnetic 
field B and V at the ground-state. Here U — 2. 




ctiemical potential: u 




FIG. 5: The fermionic concurrence as a function of chemical 
potential fi (u instead of /i as an indictor for axes) and V at 
r = 0.1. Here U = 2. 



FIG. 7: The fermionic concurrence as a function of magnetic 
field B and V at T = 0.1. Here (7 = 2. 



to add site energy term to the Hamiltonian jnj 



(14) 



For simplicity, we set fj, = fii — — /i2j since the main 
physics is not changed under this condition. 

We show the fermionic concurrence as a function of fi 
and V as well as its corresponding contour map in the 
ground-state for a specific value of [/ = 2 in Fig. 01 
Clearly the ground-state phase diagram can be obtained 
from the contour map. The area outside the arc is the 
spin singlet state described in Eq. ((S)), while the inside 
area is one-particle state. Thus the arc in the contour- 
map is just the level crossing point for different states. 
The level-crossing also leads to a jump on the surface of 



the concurrence in the ii—V plane. Since the experiments 
are always done at finite temperatures, it is useful to 
investigate the effect of thermal fluctuation. We show the 
dependence of the fermionic concurrence on fi and V at 
T = 0.1 in Fig. |S1 where we observe that the jump at the 
level-crossing point in Fig. 0]is eliminated by the thermal 
fluctuation. We find, at some points, eg., L'^ = 2 and 
V — 1, the concurrence is small when /i = 0. If we change 
the value of /i the thermal concurrence can be increased 
substantially. Such a role played by the local potential 
is important, because the change of /i in principle can be 
realized by the electric field so it provides a possible way 
to control the magnitude of the concurrence. 

The second case we consider is to add a Zeeman-like 



term to the Hamiltonian © 

Hzee = ~5i(nij - ni,4)/2 - BzKj - n2,i)/2. (15) 

For simplicity, we again set B = Bi = — i?2- We show 
the concurrence as a function of B and V for a specific 
value of [/ = 2 in the ground state in Fig. |51 and at 
T = 0.1 in Fig. respectively. Similar to the first case, 
the magnetic field may also induce the level-crossing in 
the ground-state, which may lead to a quantum phase 
transition. The induced level-crossing results in a jump 
on the surface of the concurrence in the B — V plane. 
Such a jump is blurred by the thermal fluctuation, as 
shown in Fig. [7| Here we have an interesting phenom- 
ena: although the fermionic concurrence is defined in a 
single sector, such as for spin-up electrons only, it can 
still be modulated by the other sector (spin-down elec- 
tron) through magnetic field. For the present case, the 
magnetic field does not always decrease the concurrence, 
it may actually increase the concurrence in some regions. 
Therefore, a nonuniform magnetic field can be in favor 
of the entanglement of formation between two qubits. 

V. SUMMARY AND ACKNOWLEDGEMENT 



It was interesting to observe that the classification of two 
different regions, one dominated by the charge-density 
correlation [| tiiO) -I- |0, ti)] while the other one domi- 
nated by the spin-density correlation [| Tj i) + I ii T)]? is 
closely related to the behavior of the fermionic concur- 
rence at the ground state. At finite temperatures, we 
also obtained the explicit expression of C for general U 
and V , and found that it is a decreasing function of the 
temperature, which leads to a condition of the threshold 
temperature where the concurrence vanishes. We then 
obtained the threshold temperature as a function of U 
and V numerically. Our method, as we pointed out in the 
introduction, could be easily extended to other electronic 
systems. Results obtained in this paper help us gain more 
insight into physical properties of strongly correlated sys- 
tems. Moreover, we studied the role of nonuniform local 
potential as well as magnetic field for the extended Hub- 
bard dimer. We found that both of them can be used to 
modulate the concurrence. Though our analysis is based 
on a specific value oiU — 2, the main physics is the same 
for other positive U. Finally, we want to point out that 
although results obtained in this paper are for a 2-site 
(dimer) system, their qualitative features should be the 
same for big size system. 



In summary, we first presented a possible way to study 
the concurrence as a measure of pairwise entanglement 
for fermionic systems. Then we studied its behavior in 
the Hubbard dimer both at zero and finite temperatures. 
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